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DESCRIPTION 



METHOD 



Inverse Interpolation, A Real Root of f (x) « (dqi or SADOl) 
Closed 

33 

S3 s address of auxiliary routine to^ compute f (x) 
10, 11, 12, 15, 1^, 15, 16 
2 or better 

[2(t+l) + n(t+6)] us, where a = number of Iterative linear 
Interpolations to find root, and t = duration in ms of 
auxiliary routine which computes f (x). 
Enter with the link address in Q and two values x. and 
x p of the argument of f (x) in memory positions 10 and 12. 
(These two values of x and the contents of A upon entry 
will be lost by the routine so if needed must be stored). 
If there is an odd number (e.g«, 1) of roots between x. 
and x p , the routine will be left with one of them in A, 
if there is an even number (e.g., 0) of roots between 
x. and x , the routine will be left with -1 = |-l| in A 
Mttltiple roots are treated as single roots «, The auxiliary 
subroutine whose address is to be entered in S3 before 
read- in must replace x in A by f (x) in A. This inverse 
interpolation routine uses storage locations 10-16 for 
quantities which must be preserved during the operation 
of the auxiliary subroutine so that other library sub- 
routines which in general use lower memory positions as 
temporary storage may be utilized in the computation of f (x* 
If (x,) and (x ) are two arguments of f(x) such 
that rtC^) ]Xf[(xg) ] < , linear interpolation is 
used to find the next approximation to a root x .. . To 
make the iterative process second order, f [ (x. ) ] x 

n 

f [x , ] is tested for sign. If it is positive, replace 



-2- 



and fICOJ by ft* n+1 ]- K " ls negative, replace 

('A^Vl' ^n 1 * f[( Vl )] ' ««*L>» ] * y (l/2) 
f [ fx ) 3 » a ' 7 *^ leave (x- ) unaltered* The process is 

repeated until either |(l/2) *(* n+1 )| - ° or 
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CRIER : 

eo k(h) 

S5 F 

k2 30L 
^2 3H. 
L5 10F 


'Lj 
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IOTB3 PAGE 1 i 
- Plant link addresses 






3 


50 2L 
26 S3 
10 IF 




- l/2f[(x^ ] toll 






k 


40 11F 
00 IF 


— 


Waste 






5 


L5 12F 










6 


50 5L 
26 S3 
10 IF 




- 1/2 f[(x 2 ) Q ] to 13 
1 






7 


kO 13F 


J 






8 


50 13F 
75 11F 
36 311* 




.- Test ftC*,^] rt&^o 1 






9 


L5 11F 
LO 13F 


from 29 








10 


kO l^F 
50 10F 










H- 


75 13F 
kO 16F 




— Linear interpolation 






12 


SI F 
50 12F 










13 


Jk 11F 
LO l6F 










U 


66 l^F 
S5 F 




— 






15 


1*0 l^F 




" 






16 


50 15L 
26 S3 
10 IF 




- 1/2 f (x n+1 ) to 5 






17 


kO 15F 
L3 15F 




J 





H 



LOCATIOH 


ORDER 




MOTES PAGE 2 


18 


36 30L 




J 




L5 10F 








19 


LO 12F 
M-0 16? 








20 


51 7L 
00 2F 






.-Tart 2" 58 - \(x 1 ) ji - (x 2 ) fi | 


21 


L2 16F 










36 30L 


* 






22 


50 15F 
75 11F 






-Test f(* n+1 ) fK^^l 


23 


32 25L 
L5 11F 




— 




2k 


10 IF 
kO 11F 






_ Replace l/2 f [ (j^ ) n ] by l/2 f [ (^ )J 


25 


22 27L 






26 


L5 12F 
^4-0 10F 


from 23 


Replace (x x ) n by (Xg) n 


27 


L5 13F 
i+0 11F 




Replace l/2 fKx^l by l/2 flO^J 


28 


L5 l^F 
^0 12F 


from 25 


Replace (x^ by x^ 


29 


L5 15F 
kO 13F 




Replace 1/2 ft&^y ** l/2 ^(x n+1 ) 




26 9L 




Repeat cycle 


30 


L5 ikW 


from 18, 21 


Put x . InA and leave 
n+1 




22 ( )F 


by 1 




31 


L5 32L 




Put »1 in A and leave 




22 ( )F 


by 1 




32 


80 F 


, 


1 




00 F 
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